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Chapter 1 



The Code 




Here we introduce Fisher4Cast, a code developed for general Fisher Matrix analysis. In addition to 
this manual there is also a release paper for Fisher4Cast which provides more background detail and 
cosmological applications as well as sample code for generating the plots in the release paper [1]. This 
Users' Manual explains how to install and run the code, both in the Command-line and the Graphical 
User Interface (GUI), which is coded for cosmology. The structure of the code is described in detail, as 
well as the various tests performed during the development of the code. 



1.1 Introduction 

The Fisher Matrix translates errors on observable quantities measured in a survey into constraints on 
parameters of interest in the underlying model. As such, it is the elegant way of doing propagation of 
errors to the case of multiple measurements and many parameters 

In contemporary cosmology, Fisher matrices are used to forecast parameter constraints from a proposed 
survey, and can be used to optimise future surveys (see [T] for a detailed discussion on the Fisher Matrix 
formalism). Fisher4Cast was developed with the aim of providing the community with a free, standard 
and tested tool for Fisher Matrix analysis, that is both easy to use through the Graphical User Interface, 
and yet also a robust general base-code for research. The underlying modular code of Fisher4Cast is 
completely general and is not specific to cosmology although the default setup for the GUI is intended 
for cosmology. It provides parameter error forecasts for cosmological surveys providing distance, Hubble 
expansion and growth measurements in a general, curved, FLRW background. 

The code is written in Matlab and can be redistributed and modified under the Berkley Software Distri- 
bution (BSD) license [3] which is now the required license for the Matlab File Exchange^, the primary 
repository for Fisher4Cast. Both the command line and GUI versions of the code produce plots that 
are generated directly from the program and can be easily edited in the GUI and saved in a variety of 
standard formats ('.eps', '.jpg', '.pdf etc..) for inclusion in research publications. There is also a report 
extension facility which generates either a text or I^Tj^X report complete with syntax for incorporating 

*http: / /www. mathworks.com/matlabccntral/filccxchange/ 
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the input and output in tables, matrices and figures. This means that both input and output data and 
results are easily portable from Fisher4Cast into a research publication. 

The simple start-up procedure and ease of use of the Fisher4Cast suite make it well-suited to both 
teaching and research. The input to the Fisher4Cast code can easily be changed and adapted, hence is 
can be run in large loops to explore parameter spaces, and for visualisation of the Fisher Matrix. We 
now describe the Fisher Matrix framework, outline the start-up procedure of Fisher4Cast, and describe 
the various functions and routines. A shortened version of the start-up procedure of Fisher4Cast is 
found in the Quickstart.pdf guide - which is included both as an appendix in [1] and in the bundle of 
Fishcr4Cast software. 



1.2 Getting Started 

Currently the code is available for download at one of the following websites [H [S]. Save this '.zip' file 
into the directory you want to run the Fisher4Cast suite from. 

1.2.1 The Graphical User Interface 

The GUI can be started from the Matlab editor. The file FM_GUI.m must be opened from the directory, 
and once the file is opened (click on the file icon from within the Command-line interface to open it with 
an editor) press 'F5' to run the code. This will open up the GUI screen. 

You can also launch the GUI from the command line by typing: 
»FM_GUI 

The output data will not be saved into the workspace, but the 'Saving Features' button allows one to 
save the input and output from any particular run in text or I^T^^JC code. 

The Basic Layout Explained 

We describe the basic layout of the GUI, and illustrate the various actions with screenshots taken of a 
working GUI. 

The GUI has three main sections. The section on the top left controls the input to the GUI. The bottom 
left panel controls the things one might like to use in the analysis and the parameters you are interested 
in plotting. In Figure (jl.ip we show the initial GUI screen, highlighting the observable about to be 
used (here the Growth function G{z)), and the cosmological parameters relevant to the analysis (which 
are the wq and Wa coefficients in the Chevallier-Polarski-Linder paramcterisation of dark energy [12i ?] 
- see Eq. (|1.3p ). The specific cosmological example is described in detail in [Tj, which contains the 
set of analytical derivatives used in Fisher4Cast. The right-hand side of the GUI controls the plotting 
commands for the ellipse. The various actions used to control the output are described below. 

Changing the Input Structure 

In order to compute Fisher ellipses for different input structures, one can cither choose from a drop down 
list of default example structures contained within the distribution (as shown in Figure (jl.2p ') or one can 
generate a unique input structure. This file can then be loaded to the GUI, which must be given as '.m' 
file. You can obviously also just edit the input parameters in the GUI after the default input has been 
loaded or alternately you can edit the input file (eg Cooray_et_al_2004.m). 
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Figure 1.1: The start-up screen of the Fisher4Cast GUI. 
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Figure 1.2: Changing the default input structure from the drop-down menu. 



Floating Help 



Floating help is provided with the Fisher4Cast GUI for most commands. The floating help is activated 
by moving the mouse pointer over the button or parameter on the GUI and leaving it there for a few 
seconds. This generates a screen prompt, which pops up and gives information about the function of the 
button or parameter in question. Figure (|1.3|) shows this help prompt for the 'Run' button on the GUI. 




Figure 1.3: The floating help for the 'Run' button. 



Running risher4Cast 



Once satisfied with the observables considered and the parameters of interest, pressing the 'Run' button 
will execute the code. A box will pop up that will state that the code is running, and an error ellipse will 
appear when the code has finished running. This is shown in Figure (|1.4p . 




Figure 1.4: Running the Fishcr4Cast from the GUI. 



Errors in the Input 



When the 'Run' button is pushed, the GUI first calls the FM_errorchecker.m function with the input 
supplied. This checks for the input files, checks that the data vector (e.g. redshifts at which one has 
measurements of the Hubble parameter) and error vectors (e.g. the fractional errors on the Hubble 
parameter, an / H, at the redshifts above) are the same length and performs other consistency checks. 
Should any of these tests fail, an error box will appear explaining which errors to fix before calling the 
GUI again. A log file of these errors is created in the same directory the GUI is being run in, and is 
called 'log. mat'. Loading and reading this log file is described in Section [1.2.31 Figure l|1.5p shows the 
error dialogue box indicating that a single error has been found. 
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Figure 1.5: The GUI error message box. The errors are detailed in 'log. mat'. 



The Fisher Ellipse 



Once the code is running smoothly, the resulting Fisher error ellipse is plotted. This is shown in Fig- 
ure (HI]). 





Figure 1.6: The Fisher ellipse resulting from a run of Fisher4Cast. 



Plotting more than one Ellipse 



Should one want to superimpose more than one ellipse, click the 'Hold on' button. This works both for 
the line and the area (although the same line and area fill properties will be used for both ellipses - see 
the below item for discussion of changing the colour of the area fill). Figure (|1.7p shows the resulting 
ellipse for two observables, G{z),dA{z)- 




Figure 1.7: Multiple ellipses arc plotted in one figure. 



Area fill 



Clicking on this button yields a filled error ellipse. Once it is clicked a colour must be selected from the 
menu on the left pop-up box. Note that should more than one error ellipse be plotted later, this area fill 
box must be ticked and un-ticked again to change the colour, otherwise the same colour will be used for 
all filled ellipses. This box is shown in Figure ()1.8p . 




Figure 1.8: The Area Fill option with colour selection. 



Importing Data 



The input data can also be imported from a file - either as the redshift vector (the data), the error vector 
or the matrix of prior information on the cosmological parameters. This can be done by clicking the 
relevant 'Browse' buttons on the GUI. This brings up a screen in which one can either load the data from 
file, or from the clipboard, in which case the data is cut and paste into the GUI fields. Figure (|1.9[) shows 
the screens for the loading of data from a file in the directory. In addition there is a check-box which 
specifies whether or not to use the prior matrix. 
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Figure 1.9: Loading data into the GUI. 



Multiple a 

It is possible to plot the ellipses for multiple confidence levels (i.e. 67%, 95%, 99% specified by 1—, 2—, and 3- 
a respectively). This is done via a drop-down menu on the right-hand side of the GUI, and is illustrated 
in Figure (jl.lOp . It is worth noting that the 'Hold on' was used to generate this plot. 
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Figure 1.10: Ellipses with many a. 



Different Definitions of the Figure of Merit 



Fisher4Cast provides various Figures of Merit in order to compare different surveys. These are defined 
in detail in Section fl. 4. 41 Figure (|1.15[) illustrates how these various FoMs are accessible in the GUI. 




Figure 1.11: Various Figures of Merit a are calculated in Fisher4Cast. 



Controlling Output 

The buttons on the right-hand side of the GUI all control the output specifications of the ellipse, such 
as the limits of the x and y axis, the line style and colour of the ellipse, and whether or not to have a 
grid on (over or under) the data. This is designed for maximum flexibility in representing the ellipses in 
a unique and distinguishable way. The axis labels can also be modified using the 'Edit Axis' button. 



Saving the Plot 

Once satisfied with the ellipses plotted, figures are saved by clicking on the 'Saving Features' menu and 
selecting the 'Save Plot' option. This will bring up a window to save the figure to a particular directory, in 
a selected file format ('.eps', '.fig', '.png', '.pdf etc.), as illustrated in Figure (|1.12p . No legend information 
(for example which colours correspond to which ellipses) will be saved. If legend entries are required, the 
figure must be saved as a '.fig' file, opened with Matlab, with the required legend entries added later. In 
addition to the simple saving of the figure, the reporting features of the Fisher4Cast code are accessed 
through this 'Saving Features' menu. In short, they provide the opportunity to save a report of the inputs 
and the resulting outputs and Fisher matrices from the survey, in either ASCII text or MJ^JC format. 
When clicking the 'Text report' feature the user is prompted to save the resulting '.txt' file. Similarly if 
one chooses the I^T[t;X report, both a '.tcx' report and an Encapsulated Postscript File will be generated, 
and the user is prompted to save both the figure and the M^TjX document. These reporting features are 
discussed in more detail in Section fl. 6. 21 
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Figure 1.12: Saving the figure of your Fislier4Cast ellipse. 



Skins 



The GUI is available in a variety of skins and backgrounds. These can be chosen from a drop-down list 
(consisting of both colour schemes and background images [6l[7l[8]); additional background images can 
be loaded by the user. This is shown in Figure p.l3p . 




Figure 1.13: Changing the skins and background images. 



Fisher4Cast Menu 



A Fisher4Cast menu is defined in the top left-hand corner of the Fisher4Cast GUI. From this drop-down 
menu one can access the Readme file of the code suite, the Users' manual and Quickstart Guide 
for easy reference, and the version history of the code. The BSD licence [3] for the Fisher4Cast suite 
is also available from the drop-down list. This list is illustrated in Figure (|1.15p . In addition to the 
Fisher4Cast menu, there is a menu with information on the extensions available for Fisher4Cast. In 
future releases of the code this menu will select the extensions themselves, at present it provides the 
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Figure 1.14: The Fisher4Cast drop-down menu with information on the code and version. 



Interactive Plotting 



Interactive 'point-and-click' plotting is available in Fisher4Cast Version 2.0, available by selecting the 
'Activate Interactive Plotting' option from the 'Fisher4Cast Extension' menu. Once selected the user 
interactively sets the values for parameters being plotted by clicking on the plotting area of the GUI. The 
arrow is activated for the first click on the plot area, the next click will run the code and produce the 
appropriate ellipse. The values selected will be displayed in the parameter input section of the GUI, the 
same as they would have been if manually entered. Care should be taken not to step to very unphysical 
values of the parameters, e.g. very positive values of wq. 
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Figure 1.15: The Interactive Plotting Feature in the Fisher4Cast code. This allows you to click on the 
parameter plane to change the assumed fiducial model. 



1.2.2 The Command Line 



Running the Code 

Open your version of Matlab and change the working directory to be the same as where you saved 
Fisher4Cast in. To run the code from the command hne with one of the standard test input structures 
supphed, type: 

»output = FM_run(Cooray_et_al_2004) 

This wiU caU the code using the pre-supphed test input data (Cooray_et_al_2004) and then generate 
an error eUipse plot for the parameters and obscrvablcs supphed in the chosen input. AU the relevant 
generated output is written to the output structure. You can see the range of outputs to access by typing: 

>>output 

and then examine each output individually by specifying it exactly. For example: 

>>output .marginalised_matrix 

will access the marginalised_matrix field in the output structure. It is worth noting that each '.' denotes 
another sub-level in the input structure. 

Example input files arc supplied as a template for generating new input files with your own customised 
parameters and values. All fields specified in the example inputs must be specified in any user-defined 
example input. These are outlined in Section [1.4. 21 

The code can also be run from the Matlab editor. Once the code is opened (open it from inside the 
Matlab window), pressing 'F5' will run the code. Note that if the code is run from the Editor it will call 
the default input structure, which is the Cooray_et_al_2004.m file. This is an example file containing 
input data from the paper by Cooray et al. [5] . This output can be directly compared to that of Figure 1 
of that paper. If your output compares correctly, you have a working installation of the code. Another 
input available is Seo_Eisenstein_2003.m [TO] . 



1.2.3 FM_err or checker 

The error-checker function acts 'behind-the-scenes' to check that the input structure and all the required 
variables are correct before executing the code. It can be run directly by using the command: 

>>FM_errorcheck(FM_initialise) 

where FMJnitialise is the specific function to initialise the input structure. The error checker validates, 
among other things, that all the derivative functions (whether analytical or numerical derivatives are 
going to be implemented) do in fact exist and that the data and corresponding variances vectors are the 
same length. This error checker is continually being updated to facilitate ease of use of the code. All 
error and checking messages are displayed to the screen and are also saved in the Matlab file 'log. mat' 
which can be loaded and examined at a later stage by invoking the following command: 



>>load(log.mat) 



1.3 Flowchart 



Throughout the code structures are used to allow different sub-parts of the general code access to the 
data. These structures are defined as global variables. The structures containing information on the 
input for the code are either defined at the beginning or loaded from file. All output structures can be 
saved for later use. We now outline the general framework of the code and describe the cosmological 
example specifically coded in Fisher4Cast. The flowchart shown in Figure (|1.16p summarises the layout 
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Figure 1.16: The fiowchart of the code showing the processes, decisions and storage of data. For a key to 
the symbols see Figure (|1.17p . 



and structure of the code, and the symbols are summarised in Figure (|1.17p . 
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Figure 1.17: This key shows the symbols contained in the flowchart, Figure ()1.16p . They are from, left to 
right, the begin and terminate indicator; a simple processing function which would generally return an 
output; an if statement or for loop; an input process function designed to be edited and changed as per 
the user specifications and lastly a stored structure for either input or output and passed globally for use 
throughout the code. 

1.4 Components of the Code 

We now discuss the various components of the code in detail. With each section we give a subsection of 
the flowchart to highlight the position in the flow of information through the code structure. 



1.4.1 FM_run 



FM_run.m is the general wrapper of the code. In order to make the code clear and easily editable, 
all main processes are called from this general function, and it is where all data storage occurs. Links 
to separate functions for specific calculations are documented in the code. From the command line this 
code is called with one argument, namely the specific function that initialises the input structure. This 
is unique to each example. As outlined in the section on implementing the code, if no argument is given 
the code is run with a pre-defined function, given by Cooray_et_al_2004.m, which gives the parameters 
for a redshift survey as outlined by Cooray et al. [9] . 



1.4.2 FMJnitialise 



FM initialise 



>( FM_input 



This function initialises the input used throughout the processing of the Fisher4Cast code. The values, 
names and areas of interest are specified here. It is called by FM_run.m to set the initial values for 
the input structure, which is then passed globally to all other parts of the code. Examples of the 
default initialising functions provided are Cooray _et_al_2004.m, Seo_Eisenstein_2003.m [9l[Tn|. It 
is important to note that the format of the initialising function must be kept constant - in other words 
the same input must be specified in each initialising function - the code expects values for certain entries 
in the input structure. The entries are as follows: 



• input. function_names - A cell of strings containing the specific filenames of the analytical deriva- 
tives. Note in the coded cosmological example that no analytical derivative function is specified for 
the growth function, derivatives are only taken numerically. 



• input. observable_names - A cell of strings with the names of the obscrvables. 

• input. observableJndex - A vector of the indices corresponding to the observable_names you are 
interested in, eg [2 3] would imply that your considering the second and third observables listed in 
observable_names . 

• input. dataji} - These are the row vectors of the data for each of the respective observables (indexed 
again from beginning to end by i). In the cosmological example input. datajl} would be the 
rcdshifts at which you have measurements of the Hubble parameter, for example. 

• input. parameter _names - A cell of strings containing the names of the parameters you can include 
for consideration to generate Fisher Ellipses. 

• input. base_parameters - A row vector of the parameter values (they must be specified with the 
same order as the parameter _names vector). This is the model assumed to be true in the analysis, 
the Fisher Matrix is taken around this fiducial model. 

• input. prior_matrix - The prior matrix for the parameters taken from previous surveys etc. The 
order of the matrix columns and rows correspond to the respective parameters listed as your pa- 
rameter _names . 

• input. parameters_to_plot - A row vector of the indices of the specific parameters you want to 
plot. If one index is given then a likelihood function for that parameter is plotted and if two are 
specified then an error ellipses is plotted. Selecting more than two parameters will produce an error 
message, as Fisher4Cast is only coded for up to 2-dimensional error contours. 

• input. num_parameters - This is a derived value and is given by the length of the number of 
parameters you are considering in parameters_to_plot. 

• input. num_observables - This is a derived value and is generated from the number of observables 
under consideration in observableJndex. 

• input. error{i} - The fractional error on the data from the observables (crx° /X"). It is key that 
there are as many error entries as there are observables you are considering (i.e. input. errorjl} 
gives the error on your measurements of the Hubble parameter, measured at input. data{l}). The 
entry can either be a row vector, in the case of uncorrelated observables (this vector is converted 
to a diagonal matrix in the code) or a covariance matrix. 

• input. numderiv. flag - A logical entry is expected here for each observable, should you wish to 
use numerical derivatives. 

• input. numderiv. f - Single string entries which arc combined into a struct later. These entries are 
only required if you have specified that you would like numerical derivatives for your observables. 
They give the function name of the function (say g.m) of which you arc taking derivatives. 



1.4.3 The Derivative Loop 

The code now runs various operations in a loop over the specific observables. For each observable it checks 
whether numerical or analytical derivatives are to be used by checking the numerical fiag as specified in 
the input structure (see the above discussion on the numderiv. flag for the input structure). Both the 
analytical and numerical derivatives return a matrix of derivatives for all the parameters and obscrvables 
as well as a vector of the function evaluated at the data points specified. The specific details of the 
numerical and analytical derivative codes are outlined in the proceeding discussion. Once the selected 
derivative process is completed, the relevant output is stored in the output structure. 
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Numerical Derivatives 



The numerical derivative code will calculate the numerical derivatives of any function (say g.m) provided 
that the function is specified as a function of the input parameters (i.e. g = g(d, 9x, 0b, ■■■)), by calling on 
FM_process_numeric.m which in turn calls FM_num_deriv.m and passes it the name of the function 
you wish to take derivatives of. 

The standard numerical derivative algorithm used is known as the complex-step method [llj : 

g{d,9A + ih,eB,...y 
h \ ' 

where Im represents the imaginary part of the argument, and = — 1 as usual. This method is a second 
order accurate formula and is not subject to subtractive cancellation. Unlike the finite-difference method 
an arbitrarily small step-size can be chosen and therefore the complex-step method can achieve near 
analytical accuracy. 

In addition, the simple double-sided central derivative is coded in the FM_num_deriv.m function. In 
order to use this algorithm the user must change the method field inside the derivative function from 
'complex' to 'central'. In this case the gradient is then calculated as 

dg ^ g(d, + h, 0b, ...) - g(d, 9a - h, 0b, ...) 
d0A 2h 

This is then iterated until the gradient converges for the parameter. Note that the convergence criterion 
is quite stringent - and an error message will result if there are possible convergence issues. However this 
criterion can be relaxed by changing the settings in the FM_num_deriv.m code. 

Once the derivatives are saved the Fisher Matrix must be calculated for this observable. This is done 
by first calculating the data covariance matrix for the observable. This is done in FM_covariance.m 
which is passed the function value and the index of the observable. The code checks if the error entry is 
a covariance matrix (in the case of correlated observables) or a vector in the uncorrclatcd case. It then 
calculates the covariance matrix by multiplying the variance with the function value at the data points 
considered. 



dg 
OOa 



= Im 



FM_matrix.m then produces a Fisher Matrix (F) from the covariance matrix (C) and the derivative 
matrix (V) using matrix multiphcation as F = V'^C~^V. 



Analytical Derivatives 



FM_analytic_deriv_$name 
1 

FM_process_analytic 



The analytical derivatives are specific to each user and example. If one knows the analytical form of both 
the function and of the Fisher derivatives, one can include these functions explicitly. The only conditions 
on these functions are that they must be of the form g = c?(d, 6) and must return as output a matrix of 
Fisher derivatives dg/dO and a vector of the function itself evaluated at the data points d given in the 
input. data{i}. These derivative functions are supplied for the Hubble parameter and angular diameter 
distance as FM_analytic_deriv_l.m and FM_analytic_deriv_2.m respectively. The Fisher derivatives 
of the angular diameter distance with respect to the cosmological parameter 17^ must be taken as Taylor 
series expansion when 51^ — > (see [T] for the full set of derivatives in Fisher4Cast). 

As in the numerical derivative case, once the derivatives are saved the Fisher Matrix must be calcu- 
lated for this observable. This is done by calculating the data covariance matrix for the observable in 
FM_covariance.m which is passed the function value and the index a of the observable. The code 
checks if the errors are specified a covariance matrix (in the case of correlated observables) or a vector 
in the uncorrelated case. It then calculates the covariance matrix by multiplying the variance with the 
function value at the data points considered. 

FM_matrix.m then produces a Fisher Matrix [F) from the covariance matrix (C) and the derivative 
matrix (V) using matrix multiplication as F = V'^C~^V . 



1.4.4 Final processing 
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F]V[_sum.m collates all the derivative matrices from the previous steps and sums them together to form 
a full Fisher Matrix. The individual Fisher matrices for each observable are added to the prior matrix (as 
specified in the input structure). This complete Fisher Matrix is general and is assigned to the output 
structure for future reference. 

FM_marginalise.m produces a marginalised Fisher Matrix (say F). It takes the parameters you are 
interested in (specified as parameters_to_plot in the input structure) and shuffles the Fisher Matrix 
into a block form. It then performs matrix multiplication on the blocks to produced the marginalised 
Fisher Matrix, which is also assigned to the output structure. 

F]V[_output_fom then produces the appropriate error for the likelihood case and a range of FoMs, listed 
below, for the case of an ellipse (e.g. if the length of parameters_to_plot in the input structure is 
two then an ellipse will be plotted). These Figures of Merit arc then stored in the output structure. 
Fisher4Cast includes the standard FoMs as well as some new ones available through the GUI and command 
line. Although some of the Figures of Merit are only defined for the error ellipse in the wo~Wa plane, where 
Wo, Wa are the coefficients in the CPL [121 [13] paramcterisation of the equation of state of dark energy (see 



for e.g. [11]), the FoMs in Fisher4Cast are calculated by the code for the pair of cosmological parameters 
being considered rather than the full 5-D matrix. We briefly outline the FoMs used in Fisher4Cast: 

• DETF 

This Figure of Merit in the Report of the Dark Energy Task Force [HI is defined to be the reciprocal 
of the area of the 2 — cr error ellipse in the wq — Wa plane of the CPL dark energy parameterisation 
[12l[13]. This is, equal to det{F'^^'^)/{TrV6J7). Unfortunately the DETF report does not appear to 
use this definition, and instead quotes det(i^^/^), which is the inverse of the 1—a ellipse in units of the 
area of the unit circle. Because of the benefits of the geometric interpretation Fishcr4Cast returns 
the true inverse area of the 2 — cr ellipse. To convert from one DETF FoM to the other, one should 
multiply the Fisher4Cast DETF output by nVKU ~ 7.8 

• Area^j':^. 

This Figure of Merit is the reciprocal of the 1 — a error ellipse area in the parameter plane currently 
plotted, i.e. det{F^/^)/{TrV2M) 

• Areai_CT 

Simply the inverse of the previous FoM. 

• TrC 

This FoM is defined as the trace of the covariance matrix of the data, C = F^-*-, estimated as 
the inverse of the marginalised Fisher Matrix. This FoM is simply the sum of the squares of the 
marginalised errors on each parameter. 

• Sab ^ab 

This FoM is defined as the sum of the squares of the entries of the whole covariance matrix, C = F^-*- . 
Unlike the previous definition this FoM is sensitive to the off-diagonal components of the covariance 
matrix as well as the diagonal components. 
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The data from all parts of the code are saved in the output structure. The structure formalism in 
Matlab means that each '.' indicates a further sub-level in the structure. Entries in the structure 
are of mixed type (i.e. output. function_value is a cell of vectors, one for each observable, while 
output. function_derivative is a cell of matrices of derivatives, again with one matrix of derivatives 
for each observable). By the end of the execution of FM_run.m the output structure should have the 
following entries: 



• output. function_value - This is a cell which contains a set of vectors for each of the observables 
considered. If in the specific run of the code you have only calculated an error ellipse for say one 
out of three observables then the rest of the entries are empty vectors. 

• output. function_derivative - This cell now contains matrices of the Fisher derivatives for the 
observable. Again, the entries of observables you are not considering will result in empty matrices. 

• output. data_covariance - This cell contains the calculated data covariance matrix corresponding 
to each of the observables considered. 

• output. matrix - This cell contains a separate Fisher Matrix for each of the observables considered. 

• output. summedjnatrix - This cell contains the sum of the Fisher matrices for each observable, 
and the prior information matrix, if included. 



• output. marginalised_matrix - The marginalised Fisher Matrix given here depends on which 
parameters are of interest in each run of the code. The marginahsation via matrix multiphcation 
is outhned in Section fl. 4. 41 

• output. fom - This vector contains either a single entry (1 — cr error), in the case where a one- 
dimensional likelihood function a parameter dp^ (for example) is being considered or an array of 
different FoMs when an ellipse of two parameters is being plotted. These are each explained in 
above in 11.4.41 



1.4.5 Generating plots 
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FM_generate_plot calls the either FM_plot_ellipse.m or FM_plot_likelihood.m depending on whether 
a 1— D likelihood or an ellipse is required (whether one or two parameters are specified in the parame- 
ters_of_interest field in the input structure. The style of the plot is controlled by the FM_plot_specifications.m 
file, which controls variables such as the line style, the colour of the lines, the resolution of the grid and 
the contour level (for example 1 — cr, 2 — cr). Similarly the file FM_£Lx;is_specifications.m controls the 
X and y labels and the range of the plot that will be generated. 

Lastly FM_save_struct is called to save the input and output structures with a user specified filename. 
One could invoke this function from the command line: 

>>FM_save_struct ( ' saved_f ilename ' , input , output) 

where input and output correspond to the structures that are being saved as saved_filename-01- 
Nov-2008.mat. The date on which the structure is saved is appended to the end of the filename. If a 
filename is not specified then the default name of FM_saved_data is used. It is important to note that 
the function overwrites existing files with the same name and date and no warning is given. Care should 
thus be taken to ensure different names are specified when saving important data on the same day. The 
structures can be loaded once again by issuing the following command: 

»load( ' saved_f ilenaine-Ol-Nov-2008 ' ) 

this will make the previously used input and output structures available in the current session. 
To end with we provide a global view of the structure of the code in Figure (|1.18p . 



1.5 Cosmological Application 

1.5.1 The Background Observables 

- Hubble parameter 

The expansion history of the Universe is described by the Hubble parameter, which is defined as 

H'iz) = HlE\z) = Hi (r!„,(l + zf + (1 - - ^lk)f{z,WQ, Wa) + r!fc(l + zf) , (1.1) 
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Figure 1.18: Schematic of the relationships between the various functions, variables and structures of 
the code. The vertical bars indicate the functions, variables and structures respectively while the arrows 
illustrate how a function may produce a variable and how that variable is in turn stored in a structure 
(left to right). Alternately the arrows can also show how a variable is retrieved from a structure and 
(possibly) used in a function. The order from top to bottom shows the chronological order in which the 
functions, variables and structures are called FM_run.m. 



with the evolution of the dark energy density, f(z) given by 

/W = exp(3£i±^rfz'). (1.2) 
Assuming the Chevalher-Polarski-Linder (CPL) expansion of the dark energy equation of state [121 [T3] : 

w{z) =Wo + Wa-^— =Wo+Wa\l~ —] , (1.3) 

1 + z \ aoj 



where ao = c/(i?ovKM) is the curvature radius of the cosmos, f(z) becomes: 

fiz) = (1 + z)3(i+"'o+-») exp |-3u;,^| (1.4) 

- Angular diameter distance 

Measurements of 'standard rulers' of known intrinsic length is widely used as a probe of the cosmology 
of the Universe. The angular diameter distance relates the angular size of an object to its known length 
to obtain a measure of the distance to the object, and given by 

dA{z) = T^-^ — i=- sinh (^/ihx{z)) , (1.5) 



where 



and E{z) is as defined in Eq. (|1.1|) . 



1 + z Ho x/nr 



X(.)=f (1.6) 



- The Growth of Structure 

The growth of structure is a potentially powerful probe of dark energy [Ml [151 [161 [iTl [181 [191 [20] . Consider 
the differential equation for the evolution of perturbations in the matter density S (assuming the pressure 
and pressure perturbations of the matter are zero - p = Sp ^ 0) [^ [^ : 

S + 2H6 = AT:Gprr,6 (1.7) 

The growth function provides the temporal evolution of these density perturbations, i.e. (5(x, z) oc G{z). 
In Fisher4Cast his is solved in a general FLRW universe, and hence there is in general no analytical 
solution to Eq. (|1.7p . Under the assumption of a flat universe and a cosmological constant (or pure 
curvature) however, the growing mode satisfies the following integral form [^ I25j : 

, _ b^lmE{z) (1 + Z')dz' 

^^'>- 2 A E{z'Y ' 

where the 5/2 coefficient is chosen to ensure that G{z) l/(l + z) as z oo This expression should 
not be used however, to compute the Fisher derivatives dG/dilk,dG/dwo or dG/dwa since all of these 
derivatives violate the validity of the equation. Instead, the growth derivatives should be computed 
numerically from the solution of the full differential equation for S(x). Rewriting the Raychaudhuri 
equation in terms of the Friedmann equation and the curvature density allows one to find an equation 
explicitly showing the curvature and dynamical dark energy contributions to the friction term: 

G" + + ^ - w{x)nM^)\ ^ - ^ (f7„(x)) ^ ^ 0, (1.9) 



where the new independent variable is a; = a/oo = 1/(1 + z), oq is the radius of curvature and r2fc(x) = 
—k/{a^x^H{xY') & Q.'D^{x) = pde{x) / Pcrit{x) are the fractions of the critical density in curvature and 
dark energy respectively. Alternatively, this can be written as a differential equation in terms of ln(a;): 

d^G , 3 fl , n,{x) \ dG 3 

' ' w(x)ilr,E[x) -J- -il„,{x)G = 0, (1.10) 
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which is the equation actually solved in Fisher4Cast since it is typically more stable numerically. Appro- 
priate initial conditions for this differential equation are set deep in the matter dominated era G{zi) = 
l,dG / d\nx{zi) ~ G{zi) for Zi > lOoQ Note that as a result, the growth solutions will be unreliable if 
w{z ^ oo) = wq + Wa > (or even if it is close to zero from below) since then there will be significant or 
even dominant early dark energy. Fisher4Cast allows the user to choose the redshift where the growth is 
normalised to unity. The Fisher derivatives all satisfy dG/dOi = at the normalisation redshift. 



1.5.2 Alternative Dark Energy Parametrisations 

The Fisher4Cast GUI is hard-coded for three cosmological observables {H,dA, and G), assuming the 
Chevallier-Polarski-Linder (CPL) parameterisation [T^lIISI with parameters {wq, Wa) - see Eq. (|1.3p . This 
is true of both the functions themselves, and the analytical derivatives included in the Fisher4Cast suite. 
The general framework of Fisher4Cast, however, means that one is not restricted to this parametrisation. 
As can be seen from Eq. (jl.ip . (|1.5p and (|1.10p . the dark energy equation of state enters in the cosmological 
observables through the evolution of the dark energy, via f{z), defined in Eq. (|1.4p . Hence for any given 
one only needs to specify (in the input structure) the names of the functions that will replace the 
current versions of FM_function_l.m (H), FM_function_2.m (DA) and FM_function_3.m (G). The 
same is true for the derivatives - either they can be coded analytically for the particular parametrisation 
of dark energy, or the derivatives will be evaluated numerically from the functions specified in the input 
structure. 

As a caveat, the GUI can only be used if the new parametrisation of dark energy still contains only two 
coefficients. If this is not the case, Fisher4Cast must be run from the command line version. 



1.6 Extensions 

The general philosophy of Fisher4Cast was to make it as easy as possible to mould and extend to the 
needs of a general user. In line with this philosophy we have introduced extensions as a means to add 
functionality and customisation to the existing Fisher4Cast suite. As a design philosophy for future 
extensions we envisage that extensions do not alter the core functions of the code but rather access the 
input, output or other modular core-functions. This will enable a large community of contributors to add 
and make available their own specific extensions while ensuring that the robust design features of the 
core code remain intact. 

An important element in ensuring the success of shared extensions is that all contributors have a good 
appreciation of the structures used in Fisher4Cast while also documenting and commenting the details of 
their code thoroughly, including the purpose of the extension, the required input and output produced, 
which files or structures the extension interacts with from Fisher4Cast and whether it is run from the 
command line or GUI. 

The latest extensions included in this release are listed below. 



1.6.1 Obtaining Baryon Acoustic Oscillation Errors from Survey Parameters 

Two modules are included that calculate errors on the Hubble parameter and angular diameter distance 
in BAO surveys characterised in input structures of survey parameters. The provided codes are extensions 
to Fisher4Cast, and should be placed in the same folder as the main code suite so that they can access 
the required elements of the Fisher4Cast suite. 



tin particular one can compare the results with the growth code available at 
|http : //gy udon . as . utexas . edu/ ~komat su / CRL / [ 



The first of these extensions, EXT_FF_Blake_etal2005 uses tlie fitting formulae of Blake et al. [26] 
to calculate the errors on the Hubble parameter and angular diameter distance given certain survey 
specifications, such as the survey area and the redshifts used for the measurements of H and d^. These 
are either given as central redshift bins or as the edges of redshift bins. The galaxy number density is 
also required, and is expected in units of 10^'^Mpc~'^/i^. 

All files associated with this module have the same prefix to identify them as external modules for 
the fitting formula of Blake et. al . The fitting formulae contain coefficients specific to either 
photometric or spectroscopic surveys, and hence it must be specified which survey one is considering. 
The input parameters to this module are supplied in an input structure, which is explained in the 
EXT_FF_Blake_etal2005JReadme.txt file. They are: 

• Input_survey.base_parameters: The fiducial values of the cosmological parameters are specified 
here, as {Ho,fl„i,^K-,wo,Wa), in the same way as the input parameters are defined in the main 
Fisher4Cast code. 

• Input_survey.surv_type: This specifies either a photometric or spectroscopic survey by either setting 
it to 'spec' or 'phot'. This defines which of the sets of coefficients to use in the fitting formulae. 

• Input_survey.z_type: This field indicates that the redshift data being used are either given at the 
'edge' of the bins or as the 'central' redshifts. If Input_survey.z_type='central' then an additional 
input for width of the redshift bins, dz, is required. 

• Input_survey.area: This field specifies the area of the survey in units of 1000 square degrees. 

• Input_survey.n: This gives the number density of galaxies and is measured in units of 10~''h'^Mpc~'^ 

• Input_survey.vecH and Input_survey.vecDA: These entries give the redshift vectors of the survey for 
the Hubble parameter H{z) and angular diameter distance. dAiz) 

• Input_survey.biasH and Input_survey.biasDA : These entries specify the bias for the Hubble pa- 
rameter and angular diameter distance. These two vectors should have the same length as that of 
'Input_survey.vecH' and 'Input_survey.vecDA' respectively. 

An example of a default input structure supplied in the file EXT_FF_Blake_etal2005 Jnput.m. The 

command used to run the Blake et al. [26| fitting formula is given as: 

>> [z_H_central ,z_DA_central , vol_H, vol_DA,sigH, sigDA] = ... 
»EXT_FF_Blake_etal2005_Main(Input_survey) 

Should no input structure be specified when calling the code, the default input (EXT_FF_Blake_etal2005 Jnput.m) 
is assumed. The module comprises of various smaller modules to compute the oscillation scales in the 
radial and transverse direction, and hence the errors on the angular diameter distance and Hubble param- 
eter. Figure (jl.lOp outlines the procedures in the module. The module outputs the following parameters: 

• z_H_central: A vector of the central redshifts of the H{z) data bins, related to the Input_survey.vecH 
fields specified in the input structure. 

• z_DA_central: A vector of the central redshifts of the dA{z) data bins, related to the Input_survey.vecH 
fields specified in the input structure. 

• vol_H: A vector of the volume of the redshift bins for the Hubble parameter. 

• voLDA: A vector of the volume of the redshift bins for the angular diameter distance. 



• sigH: A vector of the fractional error on the Hubble parameter: cjh/H. For percentage errors, 
multiply this by 100. 
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Figure 1.19: The flow of inputs and functions within the Blake et al. [26] extension module to Fishcr4Cast. 



• sigH: A vector of the fractional error on the angular diameter distance: o-dA/d-A- For percentage 
errors, multiply this by 100. 

The module EXT_FF_SeoEisenstein2007 also computes the errors on the Hubble parameter and an- 
gular diameter distance using the prescription set out in [27| and the sound horizon scale as given in |28| . 
This module contains a wrapper to call the Matlab version of the C code of Sco & Eiscnstein [^S] . In this 
module the code docs not need to be in the same directory as the Fisher4Cast suite, and runs completely 
independently of Fishcr4Cast. 

The module also takes an input structure (a default input structure with all fields specified is given in 
EXT_FF_SeoEisenstein2007Jnput.m) with the following parameters defined: 

• Input_survey.number_density: The galaxy number density in units of h'^ Mpc^'^ 

• nput_survey.wmap: A flag (1 or 3) specifying which WMAP power spectrum to use. 

• Input_survey.sigma8: The value of erg, the amplitude of linear clustering on a scale of 8 Mpc. 

• Input_survey.Sigma_z The line of sight real mean squared comoving distance error due to redshift 
uncertainties. 

• Input_survey.beta: The redshift distortion parameter is entered. 

• Input_survey. volume: The survey volume in units of h^"^ Gpc'^. 

Sigma_perp (the transverse rms Lagrangian displacement) and Sigma_par (the radial displacement) are 
calculated and saved to the input structure. The module is comprised of smaller functions; the flowchart 
of the module is shown in Figure p.20p . 

As in the case of the Blake et al. module, the Seo and Eisenstein module code can be called from the 
command line using: 

>> [Drms , Hrms , r , Rrms] = EXT_FF_SeoEisenstein2007_Main (Input_survey) 

The outputs of the code are the root mean square error on D/ s and i/s, where s is the oscillation scale. 
These are both given as fractional error, for percentage error multiplied by 100. In addition the correlation 
coefficient between D and H is given (r) , and the diagonal entry in the covariance matrix between D and 
H. 

To be sure that the extensions have access to the functions contained in Fisher4Cast we need to be 
sure that the extensions are either placed directly in the same folder or the path is specified to both 
Fisher4Cast and the extensions. One can use the path command to do this: 

»path (path , ' /path-to-f older/Fisher4Cast-v2 . ' ) 
»path(path, ' /path-to-f older/EXT_FF_Blake_etal2005 ' ) 

where the 'path-to-folder' is the path specifying the directory where the extension or Fisher4Cast code 
is kept on your local computer. Similarly, to run the Seo and Eisenstein |27| module, you will need to 
ensure that the respective extension is either in the same directory or the path is specified: 

»path (path , ' /path-to-f older /Fisher4Cast-v2 . ' ) 
»path(path, ' /path-to-f older/EXT_FF_SeoEisenstein2007' ) 
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Figure 1.20: Flowchart of the code in the Seo and Eisenstein BAO extension module to Fisher4Cast. 



1.6.2 Reporting Features for the Fisher4Cast Suite 

Two extension modules have been included to provide reports of the input and output structures during 
a run of Fisher4Cast. These reports can either generate an ASCII text file (.txt) or a M^]X file (.tex) 
which detail all the input and output produced by Fisher4Cast. 

In the case of the I^TgK reporting function the resulting .tex file can be compiled using IM^]X to produce a 
Postscript file (.ps) or Portable Document Format file (.pdf). This allows for a more polished presentation 
of the results generated from Fisher4Cast. It also includes a figure of the ellipse or likelihood plot which 
is embedded in the document. The additional benefit to generating a document in .tex format is that 
one can cut-and-pastc the I^TgK formatted syntax of the figure or any of the tabulated data for easy 
inclusion of the figure in an article or document containing the results from a run of Fisher4Cast. These 
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Figure 1.21: The figures on the left show excerpts from the .tex files that was generated from the reporting 
code while the figures on the right show the extracts of the corresponding .pdf file that was generated 
from the .tex file. 

reporting features arc accessible through the Graphical User Interface, by clicking the scroll-down menu 
bar labelled 'Saving Features' (see Figure (|1.12p ) for a scrcenshot of the drop-down menu). This opens a 
dialog box which in the case of the text report prompts the user for the .txt filename that it should be 



saved as. Upon choosing a M^r^X report, two dialogue boxes are opened and the user is prompted for the 
names of both the .tex file and the .eps file. 

These two extensions can just as easily be called from the command line. To generate a text report one uses 
the functions FM_report_text.m. The user is required to supply at least an input structure to generate 
a report. This input structure can either be a default input structures, eg Cooray_et_al_2004.m, or 
a user customised input. The function FM_report_text.m then calls FM_run(input) with the same 
supplied input which then produces the relevant output structure. Both the input and output used and 
generated from Fisher4Cast are then recorded in the report. The user can also specify a filename to 
saved the report as (if no .txt extension is supplied one will be added automatically). A default name 
of 'Fisher4Cast_Report-Day-Month-Year.txt' will be used, should no name for the report be specified, 
where the Day-Month- Year are the date on which the report was generated. For example the command: 

>>FM_report_text (input , ' report_name ' ) 

will generate a report with the name, 'report_name.txt', as described above. If the same report_name is 
used, the previous report will be overwritten without warning. Please specify a unique report_name to 
ensure the report is correctly saved. 

Finally there is an option of including a specific output structure in the report function. This is useful 
when generating the report from the GUI, but care should be taken when using this option in the command 
line, as one runs the risk of generating a report where the input and output are not appropriately related. 
In other words, 

>>FM_report_text (input , ' report_iiame ' , output) 

generates a report as before with the name, 'report_name.txt', using the input supplied and assuming 
that the given output is associated with the respective input. 

Much the same as the text report, the I^T[t;X report is called using FM_report_latex.m and requires 
at least an input structure. The filename the report is to be saved as can also be specified (either with 
or without the .tex extension). The MJr^X report includes the EPS figure generated from Fisher4Cast: 
as a default the figure will be saved with the same name as the .tex, except an .eps extension. The 
default name of 'Fisher4Cast_Report-Day-Month-Year.tex' is used, should no report name be specified 
The commands 

>>FM_report_latex ( input , ' report_name ' ) 

generate a report with the name, 'report_name.tex', and a figure with the name 'report_name.eps', where 
the names overwrite any existing files of the same name. Additionally, one can use a specific figure in the 
report with the command: 

>>FM_report_latex ( input , 'report_name' , 'use_fig') 

In this case there is of course no guarantee that the figure and the output from Fisher4Cast agree. 
As in the case of the .txt report, one can specify the output structure directly with: 

>>FM_report_latex ( input , 'report_naine' , 'use_figure' .output) 

which generates a report as before with the name, 'report_name.tex', using a figure called 'use_figure.eps' 
where the output and figure are assumed to be associated with the respective input supplied. 



1.7 How to produce small, good quality postscript images for 
inclusion in I^Tf^X documents 



Much of the output from Fisher4Cast is expected to be used in research pubhcations produced using 
M^TjX, in which case Fisher4Cast figures need to be saved in '.eps' or '.ps' format. Unfortimately the 
default Matlab '.eps' and '.ps' files produced tend to be large, often several MB in sized. Apart from 
making printing slow, this is a problem when submitting papers to online archives, like the arXi"\{§, which 
has a strict file limit. 

Many of the figures in this paper exceeded 1MB after saving from Matlab, and several exceeded 2MB. 
To solve this, one can instead save the files as bitmaps and then use a utility such as jpeg2ps to convert 
the '.jpeg' file to postscript with a much smaller file siz Nevertheless, achieving good compression 
and good quality results can be tricky to achieve with Matlab figures, and we share the steps that have 
worked well, here. 

• Save the files in '.eps' format in Matlab. 

• Load the eps file into Photoshop (or equivalent such as GIMP), which will immediately request to 
rasterize the file and hence request an image size (in cm or inches) and a resolution (in dots per 
inch, dpi). 

• Choose the size to be that which you want the figure to appear in your M^TjXdocument (e.g. 10cm). 
As with all bitmaps, choose it to be the physical size you will use in the end since rescaling bitmaps 
causes blurring and poor quality. 300dpi gives good resolution. The resulting file will be huge in 
Photoshop, but do not worry since it is only an intermediate step. 

• Now save the file as a '.jpeg' with high quality factor, 80% or better. 

• Using jpeg2ps or a similar utility, convert the '.jpeg' file to '.eps'. This will add approximately 
10% to the jpeg file size due to the postscript wrapper. You should be able to achieve a 10-20 fold 
reduction in eps to eps file size with only marginal reduction in quality. 



1.8 Tests of the Code 

Various tests were performed to check the correctness and accuracy of all components of Fisher4Cast, 
namely the integration routines in Fisher4Cast, the derivatives and their validity (especially relevant for 
the growth function) and the matrix manipulation and generation of Fisher error ellipses. 

1.8.1 Integration Tests 

The integration routines in Fisher4Cast are tested by comparing them to standard results and fitting 
formulae for the angular diameter and growth function respectively. 

1.8.2 Angular Diameter Distance 

The angular diameter distance, d^(z), defined in Eq. (|1.5[) . is the ratio of an object's physical transverse 
size to its angular size (in radians). Characteristically it does not increase indefinitely as z — > cxd, rather 

brief word of caution: if you plot many ellipses in the GUI, or with multiple colour variations, and then save the 
result, Matlab will save the entire history, making the resulting file largo. It is best to decide on what combinations you 
want and then save only that figure. 
§http:/ /xxx. lanl.gov 

^This process is discussed at jhttp: / / aps.arxiv.org/help/bitmap/index| and elsewhere. 



it inverts at z ~ 1 - thereafter more distant objects actually appear larger in angular size. In the case of 
r^A = 0, one can write Eq. (|1.5p as an analytical function of redshift and cosmic parameters as [501 121) : 



_C_ 2 [2 - (1 - Z) - (2 - Qrn) \/l + »m z\ 

Ho (1 + z) 



(1.11) 



The angular diameter distances for three particular cosmologies (the same cosmologies as found as Figure 
(2) of [3T]) are shown in Figure (|1.22p . The same axis lengths and line styles are used for easy comparison 
of the two plots, which show agreement between the Fisher4Cast algorithms and the know results. 



Dimensionless Angular diameter distance 




Figure 1.22: The angular diameter distance - plotted for Universes with (il,„,17A) = (1,0) (solid), 
(0.05, 0) (dotted) and (0.2, 0.8) (dashed hue). Our results are in agreement with a similar plot. Figure (2) 



1.8.3 Growth function 



In order to test the numerical growth function, the solution is compared to other numerical solutions in 



the literature (code is available for comparison at |http: / / gyudonTas.utexas.edu / ~ komatsu/CRL/, which 
implements the growth equation given as Eq. f76) in |32| ). to analvtical approximations for particular 
cosmologies, such as Eq. (jl.Sp . or to fitting formulae, such as that originally suggested by Carroll, Press 
& Turner [33]. This fitting formula is given by G(z) = [M] where 



g{z) = 5 
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(1.12) 



and 
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The comparison between the growth function from Fisher4Cast and the fitting formula is shown in the 
left-hand panel of Figure [1.231 together with the relative difference, (G — Gfit)/G, between the two 
methods, in the right-hand panel. 
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Figure 1.23: The Fisher4Cast growth function compared to a fitting formula. The left-hand 
panel shows the growth function implemented in Fisher4Cast (blue curves) as compared to the fitting 
formula Eq. ([112)1 [35l [34] (black dots). The models represented are ($!,„, fife) = (0.3,0.7) (dashed), 
(0.3,0) (solid) and (1,0) (dot dashed). The normalised differences between the growth function used 
in Fisher4Cast and the fitting formula for the various models are of order 10"'^, which is shown in the 
residuals of the right-hand panel. 



Degeneracy Tests 

The degeneracy direction of a Fisher ellipse (we consider WQ,Wa) is a useful diagnostic of whether or 
not the Fisher ellipses are being calculated correctly. The direction of degeneracy can be computed 
analytically for a given redshift by assuming that the specific observable X" is constant at a particular 
redshift, and to then solve for Wa as a function of wq, or by computing the likelihood over a grid of wq — Wo, 
and then taking contours of the likelihood that correspond to the same cosmology as that assumed in the 
run of Fisher4Cast. This is particularly important when considering numerical derivative routines, such 
as those needed for the growth function. 



Growth function 



As a test of the correctness of the solution (and Fisher derivatives) of Eq. (ll.lOp , the Fisher ellipse from a 
survey consisting of a single measurement of the growth function at z = 3 (as computed with Fisher4Cast) 
are shown in Figure (|1.24p . overlaid with contours of the likelihood corresponding to the growth value of 
the fiducial cosmology. The agreement of the degeneracy directions indicates that numerical computation 
of both the function and its Fisher derivatives is sound. 



Hubble parameter 



In the case of the Hubble parameter one can compute the degeneracy direction analytically: consider 
a single perfect measurement H{z) at some particular redshift z. Solving for Wa in terms of wq by 
substituting Eq. ([Ti]) into pTT]) . yields 



In 

C{z) ^ - 



ln(l - - ^k) 

= 3(l + wo+Wa)ln(l + z)-3wa-^, (1.13) 




Figure 1.24: Grid test of the numerical growth derivatives evaluated at a single redshift: The 

growth function is evaluated at a single redshift z ~ I for a range of models on a grid of —1.5 < wq < 0.5 
and —l<Wa<l- The wq — Wa degeneracy direction in an assumed ACDM cosmology is then the 
iso-growth contour corresponding to ACDM (wq = —l,Wa = 0), shown here as the orange dashed line. 
The Fisher4Cast degeneracy direction (computed assuming an arbitrary value of 1.5% error on growth) 
is shown as brown solid lines. 




Figure 1.25: Lines of degeneracy — comparing the analytical degeneracy directions for the Hubble 
parameter (long dashed lines) with the degenerate error ellipse produced using Fisher4Cast (solid lines) 
for two separate surveys of 10% measurements of one at z ~ 0.5, and another atz = 2. The analytical 
solution and that of the error ellipse coincide, which is a systematic check of the computations and 
correctness of Fisher4Cast. 



hence 

_ C(z)-3(l + wo)ln(l + z) 
3(ln(l + z)-z/(l + 2)) 



(1.14) 



This degeneracy direction is shown in Figure (|1.25p for two different values of the central redshift z along 
with the ellipses from the ellipses from Fisher4Cast. The agreement between the analytical solution 



for the degeneracy direction, and the degenerate cUipsc from the Fishcr4Cast code confirms that the 
derivatives are being calculated correctly, and the matrix operations within Fisher4Cast are sound. 



Angular Diameter Distance 



Furthermore, it is possible to obtain an analytical solution for the direction of degeneracy between the 
parameters wo,Wa from measurements of the angular diameter distance, dA{z) (Eq (jl.Sp ). The dark 
energy parameters only enter Eq (jl.Sp in x, the integral of l/E, hence we need only take derivatives 
dx/dwi in order to determine the degeneracy, via: 



-Wo 



dwo 
dx/dwQ 



Wo + ' 



(1.15) 



where Q is an arbitrary constant. Figure (jl.26p once again confirms the agreement between the degeneracy 
direction of the ellipse produced using Fishcr4Cast and the analytical solution, for two values of the central 
redshift z. 




Figure 1.26: Lines of degeneracy — similarly to Figure ()1.25p . the analytical degeneracy directions 
for the angular diameter distance (long dashed line) are superimposed on the degenerate error ellipses 
produced from Fisher4Cast (solid lines) for two surveys consisting of one measurement each of (Ia with 
a 10% error at rcdshifts z — 0.5 and z ~ 2 respectively. 



Combined Degeneracies 

Combining measurements from multiple probes is key to tightening constraints on parameters in a cos- 
mological model. Summing over different obscrvablcs X" means that one effectively "combines" the 
different degeneracy directions of each observable together. Given the analytical degeneracy directions 
discussed earlier, one can investigate (qualitatively at least) how the these degeneracy directions combine 
to form the Fisher ellipse. This is illustrated for the Hubble parameter and angular diameter distance at 
a redshift of z = 0.5. The orange dashed line shows the degeneracy direction for the angular diameter 
distance, the degeneracy direction for the Hubble parameter is shown in blue - the combined ellipse in 



the WQ — Wa plane for these two observablcs (black ellipse) is a combination of both degeneracy directions. 
Qualitatively, the magnitude of the Fisher derivatives for either observable (an indication of how sensitive 
the observable X" is to the parameter 9a) determines how much its degeneracy direction contributes to 
the final ellipse. This can be used to investigate how the degeneracy directions change with redshift. 




Figure 1.27: Contributions to Fisher ellipse The separate degeneracy directions (orange dashed line 
for (Ia and (blue dashed line for H) arc plotted along with the Fisher ellipse for wq — Wa computed at 
a redshift of z ~ 0.5, illustrating how the various observables contribute to the overall direction of the 
Fisher ellipse. One can then employ this strategy at varying redshifts to track how the contributions 
change as a function of redshift, as the observables become more or less sensitive to the dark energy 
parameters wo,Wa- 
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